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ABSTRACT 

The connection between mass loss on the red giant branch (RGB) and hor- 
izontal branch (HB) morphology in globular clusters (GCs) has long been ac- 
knowledged but the mechanisms governing mass loss remains poorly understood 
from a theoretical perspective. The present study uses synthetic HB models 
to demonstrate for the first time that a-enhancement and a simple relation be- 
tween mass loss and metallicity can explain the entire range of HB morphology 
(characterized by the HB type index) observed in old, coeval GCs. The mass loss- 
metallicity relation accounts naturally for the fact that the most metal poor GCs 
([Fe/H] < -2) have redder HBs than is typical of GCs with -2 < [Fe/H] < -1.5 
without invoking younger ages. These results may prove useful in studying the 
contribution of HB stars to integrated light via stellar population synthesis. 

Subject headings: globular clusters: general — stars: horizontal branch 



1. Introduction 



Zinnl ( 119931 ) recommended that Galactic halo GCs be divided into two groups: one for 
old, roughly coeval GCs that follow the mean trend in [Fe/H] vs. HB type in the inner halo 
(the Old Halo; OH) and one for GCs that deviate from this trend by a significant amount 
(the Younger Halo; YH). Due to the lack of very metal poor GCs in the inner halo, Zinn 
used more distant GCs to extend the line below [Fe/H] < —2, see his Figure 1. 

Zinn's relation was empirical ly motivated and diff ers from la t er stu dies — such as those 
by iLee. Demarque. fc Zinnl (11994 . hereafter LDZ) and iRey et al.l (120011 ) that are based on 
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synthetic HB isochrones — in that it turns back to the red at the lowest metallicities ([Fe/H] 
< -2). Trends derived from synthetic HB calculations wit h constant mass loss (LDZ) or mass 
loss that increases with [Fe/H] and age (IRey et al.ll200ll ) indicate that HB type (as defined 
by LDZ 1994) ri ses monoto nically to +1 as [Fe/H] decreases and remains there for all lower 
values of [Fe/H] . iLed (119911 ) showed that allowing mass loss to vary with [Fe/H], in a manner 
consistent with the observed trend in GCs, resulted in better agreement between his models 
and the pulsational properties of RR Lyrae stars in u Cen than was achieved by assuming 
constant mass loss. Aside from Lee's application to uo Cen, the concept of searching for a 
mass loss relation that matches the observed trend in old, coeval GCs has not been explored. 

The purpose of this Letter is to argue that it is possible to account for the fact that the 
most metal poor GCs have HB types less than +1 without requiring them to be significantly 
younger than more metal rich GCs if mass loss on t he RGB va ries as a function of [Fe/H] in 
a manner qualitatively similar to that proposed by ILed (119911 ) . 



2. Synthetic HB Model Construction 



Rood! (119731 ) introduced an effective means of modeling the distribution of HB stars in 
GCs by assuming the stars that arrive on the HB are drawn at random from a Gaussian mass 
distribution that is truncated on both ends by physical limitations: namely the initial mass 
and the He core mass of the precursor red giant stars. The method further assumes that 
stars arrive on the HB at a constant rate and spend roughly the same amount of time there 
regardless of mass. Under these assumptions, one constructs a synthetic HB by drawing 
mass and age at random and then interpolating in a grid of He-burning evolutionary tracks 
to obtain the relevant quantities. 



M any investigators, such a slRoodl (11973 ), LDZ (1990, 1994). ICatelan fc de Freitas Pacheco 

dl994J ), IRey et all (120011 ). and bassisi et all (120041 ). have used the synthetic HB method to 



show how various model parameters influence HB morphology and the pulsational proper- 
ties of RR Lyrae stars in GCs. LDZ (1994) generated isochrones in the [Fe/H] vs. HB 
type diagram for di fferent ages and H e abundances while assuming a fixed amount of mass 
loss and dispersion. IRey et al.l (120011 ) incorporated the Reimers mass loss formula into their 
synthetic HB analysis (such that the average mass loss increased with [Fe/H] and age) and 
found that, while the sensitivity of HB morphology increased compared to the findings of 
LDZ (1994), the general trend was unchanged. These isochrones were used to argue that GCs 
with anomalously red HBs (compared to the oldest synthetic HB isochrone) are younger than 
the OH GCs: one of strongest arguments in favor of age as the second parameter governing 
HB morphology in GCs. 
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This study em ploys the He core masses, RGB tip masses, and He-burning tracks from 
Potter et al.l (120071 ) to create synthetic HB models (each consisting of 500 stars) spanning the 
range of [Fe/H] and [a/Fe] observed in the Galactic halo. Simulations of 1000-2000 synthetic 
HBs are generated by allowing [Fe/H] to vary between —2.5 and —0.5 while [a/Fe], age, mass 
loss (AM), and dispersion (ctm) are held constant or varied in a specified manner. 

Figure [1] indicates how some of the governing parameters influence the [Fe/H] vs. HB 
type trend. Panel A compares trend lines derived from a simulation with constant mass 
loss (AM=0.15 M ), dispersion (cx M =0.02 M ), and [a/Fe]=+0.4 for ages of 9, 11, and 13 
Gyr; panel B shows the same mass loss and dispersion assumed in A but at a fixed age 
(13 Gyr) with [a/Fe]=0, +0.2, +0.4, and +0.6; panel C illustrates the effect of varying the 
mass loss (AM=0, 0.1, 0.2, 0.3 M ) when the dispersion (0.02 M ), age (13 Gyr), and 
[a/Fe] (+0.4) are held constant; fi nally, panel D comp ares synthetic H Bs derived from the 
He-burning evolutionary tracks of IDotter et al.l (120071 . soli d lines) and Cassisi et al.l (120041 . 



dashed lines) at 13 Gyr with AM=0, 0.1, and 0.2 M . The lCassisi et al. 



(2004) models were 



calculated at Y=0.23 with a scaled-solar mixture for all [Fe/H] while the lDotter et al.l (120071 . 
see their Table 4) models range from Y=0.245 (at [Fe/H] =-2.5, [a/Fe]=0) to Y=0.264 (at 
[Fe/H]=— 0.5, [a/Fe]=0) so it is not surprising that the the lDotter et al.l (120071 ) models yield 
consistently bluer HB types in most regions of the diagram. Figure [1] shows that there is 
considerable degeneracy in the [Fe/H] vs. HB type diagram among age, [a/Fe], and AM. 
Though not shown, stochastic effects are also significant because even though each synthetic 
HB consists of 500 stars (a number observed only in rich GCs) there is still a range of 
possible outcomes (AHB type~0.1) for the same input parameters. Only a few tens of HB 
stars have been observed in some GCs and this, combined with stochastic effects, can make 
the parameters governing HB morphology difficult to ascertain from model comparisons. 

Based on the relationships that were found by varying individual parameters in the syn- 
thetic HB models, additional simulations were generated in order to attempt to understand 
the observed distribution of old, roughly coeval GCs in the Galactic halo. In these simula- 
tions, [Fe/H] was allowed to vary over —2.5 < [Fe/H] < —0.5 while age was held fixed at 13 
Gyr. The level of [a/Fe] was drawn at random from a flat distribution in such a way that 
the average value of [a/Fe] decreased as [Fe/H] increased. Figure [2] shows how the fraction 
of synthetic HB models with a given value of [a/Fe] is assumed to vary with [Fe/H]. This 
method of distributing [a/Fe] is ba sed on observations of the q -elements as a function of 
[Fe/H] in halo GCs and field stars (jPritzl. Venn, fc Irwinl 120051 . see their Figure 3). Mass 
loss (AM; measured in M ) was assumed to vary with [Fe/H] as shown in Figure [3J The 
minimum value was set as a function of [Fe/H] as in equation [TJ 
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AM mi Jx ) = c- 



1 + exp(^) 



Where x=[Fe/H] and a, 6, and c are free parameters. The variation of AM with [Fe/H] 
could also have been described by a broken linear relation; the use of equation [1] was chosen 
because it provides a smooth variation with [Fe/H]. With the minimum value of AM set by 
equation [TJ the actual value for a single synthetic HB realization % was set by equation [2j 



AM, = AM mm (x t ) + 1^(0,0.03)1 (2) 

Where Rg{v, z) is a random number drawn from a Gaussian mass distribution having 
a mean value y and dispersion z. Thus the value of AM was more likely to lie near the 
minimum value but, given a large enough simulation, values at the ~3 a level can be found. 
Figure [3] shows the distribution of AM as a function of [Fe/H] using equations [1] and [2] for 
a simulation of 2000 synthetic HB models. The minimum value as defined by equation [1] is 
shown as the solid line while the dashed lines represent the minimum value plus 1, 2, and 
3a. 

The value of the dispersion in the mass distribution (<tm; also measured in M & ) was set 
in a similar fashion to that for AM (equation [2]) except that the minimum value was held 
constant at 0.02 and the dispersion in Rq was 0.01. 



3. Results 



Figure H] shows the [Fe/H] vs. HB type diagrams that are obtained when the variable 
AM min given by equation [1] with a = —2, b = 0.05, and c=0.14 (solid circles) and a constant 
AM min =0.14 (open squares) are assumed. Both simulations consist of 2000 synthetic HB 
models at 13 Gyr and both used equation [2] to determine AM and au as outlined in the 
preceding section. That synthetic HB models become redder as [a/Fe] increases at a given 
[Fe/H] may be a necessary part of the explanation for why at least some of the most metal 
poor GCs have redder HBs than can be found at [Fe/H] > —2 but [a/Fe] alone is not sufficient 
to explain this behavior. Assuming a constant mass loss and the [a/Fe] distribution shown 



in Figure [2 
(1994) and 



produces t he sa me sort of [Fe/H] vs. HB type behavior demonstrated by LDZ 
Rev et al.l (120011 ). Allowing AM to vary with [Fe/H] (see eq uations [J an d HP as 



described above produces a trend that matches the original OH trend of IZinnl (119931 ) . 



The GC catalog compiled by 



Mackey &: van den Berghl (120051 ) places GCs into sub- 
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populations along the lines of IZinnl ( 1993) but with additional categ ories that reflect discov- 
er ies made in the in tervening years. iMackey fc van den Berghl (120051 ) used the HB isochrones 
of iRey et al.l (120011 ) to assign a sub-population to each GC. Many of the lowest metallicity 
GCs are considered part of the YH despite evidence that their ages a re similar to other 



metal poor OH GCs , such as M92 and NGC2419, with bluer HBs (IVandenBerd 12000 



Salaris fc Weissl 120021 ). Given the lack of a substantial age difference between the lowest 
metallicity GCs with different HB types, the metal poor GCs NGC5053, 5466, 6426, and 
7078 (M15) have been included in the comparison presented here. Figure compares a syn- 
thetic HB simulation with a-enhancement and variable mass loss with the OH GC s (fille d 
circles) and four metal poor YH GCs (open circles) from IMackey fc van den Berghl (120051 ). 
In order to avoid complications, as much as currently possible, that arise from the HB 
morphologi es of GCs with mu ltiple stellar populations, NGC1851 (IMilone et al.ll2008l ) and 
NGC2808 JPiotto et al.ll2007! ) are not shown in Figure U 



The variable mass loss suggested in the previous sectio n reproduces t he di stribution 
of old GCs i n the halo. Indeed, if the He-burning tracks of ICassisi et al.l (12004 ) are used 
instead of the lDotter et al.l (120071) tracks, approximating q-enhanc ement by [M/H] ~ [Fe/H] 
+ Iog(0.638xl0[°/ Fe ] + 0.362) JSalaris. Chieffi. & Stranierol [l993h . variable mass loss with 
the same values of a, b, and c listed at the beginning of this section produces a nearly 
identical result to that shown in Figure [5j Based on these results pertaining to mass loss and 
HB morphology as well as the evidence that the four old, metal poor GC s with relatively red 
HBs are not especially young (I VandenBergl |2000| ; ISalaris fc Weissl 120021) . it is recom mended 
that those GCs be considered part of the OH as originally intended by IZinnl (119931 ). 



Discussion and Conclusion 



Catelanl (120051 ) provides a thorough discussion of mass loss and HB morphology in GCs. 



Catelan considers a handful of empirically motivated formulae that describe the mass loss 
rate as a function of stellar parameters such as luminosit y, radius, and s urface gravity. Each 
predicts that mass loss increases smoothly with [Fe/H] (jCatelanl 120051 . Figure 4). Some of 
the formulae are consistent with the observed difference in HB morphology between second 
parameter GC pairs , but not all, and no single formula is successful for each of the GC pairs 
considered (jCatelanll2005l . Figures 10-13). There is also a discrepancy b etween the fo r mulae 
and observational efforts to measure mass lost by red giants in GCs by lOriglia et al.l (120021 ) 
who found that mass loss is not strong ly correlated with m etallicity and only occurs near 
the tip of the RGB. On the other hand, lOriglia et al.l (120071 ) found that mass loss is evident 
at the level of HB and increases with luminosity up to the RGB tip in 47Tuc. The proposed 
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mass loss-[Fe/H] relation should be straightforward to test observationally with mass loss 
estimates for a large enough GC sample. 

These examples illustrate that mass loss remains a complex and poorly understood 
process. It is by no means guaranteed that the information provided by standard stellar 
evolution models is sufficient to calculate mass loss as a function of age and metallicity. 
The observed range of HB morphologies in OH GCs (which are assumed to be roughly 
coeval) of similar metallicity suggests that more than one factor must be involved. It seems 
quite plausible that mass loss may be influenced at the stellar level (by, metallicity and 
rotation, for example) and at the cluster level (by cluster mass or concentration, for example). 
Uncertainties in composition, coming from the lack of accurate abundance determinations 
in some lesser known GCs, or the possible presence of multiple populations with distinct 
compositions in a single GC, complicate the analysis and interpretation. 

The situat ion is further compl icated by the discovery of bi- and multi-modal HBs in 
some GCs, see Ferraro et al. ( 1998T) f or exa mples and a review of the subject. It has been 
suggested by, e.g. iD'Antona fc Caloil (120081 ) and references therein, that these (and perhaps 
all) GCs are comprised of chemically distinct populations that create the complex distri- 
bution of HB stars. Complications arise even in GCs with uni-modal HBs. The analysis 
of M3 by IValcarce &: Catelanl (120081 ) indicates that while the synthetic HB method devel- 
oped by Rood ( 19731 ). particularly the assumed Gaussian mass distribution, is 'a reasonable 
first a pproximation, [it] fa i ls to account for the detailed shape of M3's HB mass distribu- 
tion.' Valcarce & Catelan (boost ) also warn that there is a mismatch between observed and 
theoretical lifetime-luminosity relations on the HB. 

In the absence of a firm theoretical foundation for mass loss in stars, it is still possible 
to further our understanding of this subject and to apply that knowledge. One application 
of the results presented here is t o the modeling of integr ated light via stellar pop u lation 
synthesis. Investigators such as iLee. Yoon. &: Led (120001 ) and iLee. Lee, fc Gibson! (120021 ) 
have explored the dependence of integrated light properties on synthetic HB morphology 
and compared them with observations of extragalactic GCs. The results presented here may 
provide useful input for future studies along these lines. 

To summarize, synthetic HB models that include the effects of a-enhancement and 
metallicity-dependent mass loss (via equations [I] and [2]) can account for the entire distribution 
of HB morphologies observed in old, coeval GCs in the Galactic halo. Essentially the same 
result is obtained from two distinct sets of He-burning tracks, suggesting this is not a spurious 
or model-dependent result. As a corollary, it is recommended that NGC5053, 5466, 6426, 
and 7078, recently associated with the YH, be considered members of the OH since their 
ages are consistent with other OH GCs and their HB morphologies can be explained without 
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invoking substantially younger ages. 
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HB Type 



Fig. 1.— The effect of varying: (A) age at constant AM=0.15, <x M =0.02, and [a/Fe]=+0.4; 
(B) [a/Fe] at constant AM=0.15, a M =Q.Q2, and 13 Gyr. (C) AM at constant a A/ =0.02, age 
of 13 Gyr, and Ja/Fe1= +0.4; and (D) AM at constant cr A/ =0.02 usi ng the He-bu r ning t racks 
of ICassisi et all tooi dashed line; Y=0.23, [a/Fe]=0, 13 Gyr) and botter eTaD <j2007L solid 
line; Y=0.245-0.264, [a/Fe]=0, 13 Gyr). 
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[oc/Fe] = 



[o/Fe] = +0.2: 



[oc/Fe] = +0.4 




Fig. 2. — The approximate fraction of all synthetic HBs with a given value of [a/Fe] as a 
function of [Fe/H] in the synthetic HB simulations. Since the number of realizations in each 
simulation was finite, the actual fractions only approach the values shown here. 
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Fig. 3. — Distribution of AM (in M ) as a function of [Fe/H]. The points are from one 
simulation of 2000 synthetic HB models. The simulation used a = —2, b = 0.05, and 
c = 0.14 in equation [H The solid line is the minimum trend (equation [1]) and the dashed 
lines are 1, 2, and 3a above the minimum. 
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HB Type 

Fig. 4. — Comparison of synthetic HB simulations for which AM min was either fixed (open 
squares) or varied according to equation [1] (solid circles). In both cases, AM was set by 
equation [2j HB type varies monotonically with [Fe/H] in the fixed AM m , n case but turns 
back to the red in the variable AM min case below [Fe/H] < —2. 
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HB Type 



Fig. 5. — The synthetic HB simulation with variable mass loss (small filled cirlces) com- 
pared to the OH GC sample (l arge filled circles) and four YH GC (large open circles) of 



Mackey fc van den Berghl (120051 ) with omissions as described in the text. The AM-[Fe/H] 



relation given by equation [T] is capable of simultaneously reproducing the OH GCs and the 
four metal poor YH GCs without assuming any variation in age. 



